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In the present paper the resonances of vibrational modes in one-dimensional ran- 
dom Morse lattice are found and analyzed. The resonance energy exchange is ob- 
served at some values of elongation. Resonance 2cji = 0J2 is investigated in details. 
The interacting modes are inequivalent: the higher-frequency mode is much more 
stable in the excited state, i.e. its life-time is larger than the life-time of lower- 
frequency mode under the resonance conditions. Simple model of two nonlinear ly 
coupled harmonic oscillators is also considered. It allows to get analytical description 
and to investigate the kinetics and the energy exchange degree vs. such parameters 
as the resonance detuning and specific energy. The very similar behavior is found 
in the Morse and the two-oscillatory models, and an excellent agreement between 
analytical and numerical results is obtained. Analogous resonance phenomena are 
also found in the random Fermi-Pasta-Ulam lattice under contraction. 
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I. INTRODUCTION 

The primary goal of Fermi, Pasta, Ulam and Tsingou (FPUT) [l| [3] was to find an energy 
sharing in one-dimensional lattices with nonlinear interaction between neighboring particles. 
The authors expected the occurrence of the statistical behavior in a system with coupled non- 
linear oscillators due to the nonlinear interaction, leading to the energy equidistribution among 
the degrees of freedom. Initially, the long-wave vibrations were excited and "Instead of a gradual 
increase of all the higher modes, the energy is exchanged, essentially, among only a certain few. 
It is, therefore, very hard to observe the rate of "thermalization" or mixing in our problem, 
and this was the initial purpose of the calculation." 

Much effort has been undertaken to explain the FPUT results. Two approaches were devel- 
oped. The first one was to analyze the dynamics of the nonlinear lattice in the continuum limit, 
which led to the discovery of solitary waves [2j . The second approach, advanced by Izrailev and 
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Chirikov, pointed to the overlap of nonlinear resonances [3( and an existence of a stochasticity 
threshold in the FPUT system (4j. For strong nonlinearities (or/and large energies) the overlap 
of nonlinear resonances results in the dynamical chaos which destroys the FPUT recurrence 
and ensures convergence to thermal equilibrium. (For more details on the FPUT problem see 
recent reviews devoted to the 50th anniversary of the celebrated paper [1] in Focus Issue: The 
Fermi- Pasta-Ulam Problem — The First Fifty Years, ed. by D.K. Campbell, P. Rosenau, and 
G.M. Zaslavsky, Chaos 15(1) (2005)). 

In the present paper we thoroughly analyze the resonance between the lowest (k = 1) and 
the next (k = 2) vibrational modes previously discovered in the random Morse lattice under 
the elongation [5j. Here we also employ the lattice with the random interaction potential 
between neighboring particles under different deformation degrees. The lattice randomness 
makes the vibrational frequencies to be irregular, and the deformation allows to achieve the 
true resonance conditions/Random FPUT .attics (with random quadrat, terms) and with 
alternative masses j^l were investigated but no resonances were directly observed. 

The model of two nonlinearly coupled harmonic oscillators with the frequencies ratio lo\ : 
UO2 = 1 : 2 is also considered. An excellent agreement in dynamics and kinetics of the Morse 
lattice and the oscillatory model is observed. 

The random ct-FPUT lattice under compression demonstrates the similar behavior as the 
Morse lattice. The compression of the a-FPUT lattice in contrast to the elongation of the 
Morse lattice is explained by the different signs of cubic terms in an expansion of the Morse 
potential (U M = x 2 - x 3 ) and in the a-FPUT lattice U = x 2 /2 + ax 3 / 3 - 



II. RESONANCES IN RANDOM MORSE LATTICE AND TWO-OSCILLATORY 

MODEL. 



A. Random Morse lattice. 

We consider the finite-length lattice of iV particles interacting via the Morse potential with 
the rigid boundaries. The potential energy has the form: 

N+l 

Up = y^Ai U (yi)i Vi = Xi-Xi-i (xo = 0, X N+ i = L), (1) 

i=l 

where Xj is the displacement of ith particle from the equilibrium and U(y) = [1 — exp(— y)] 2 
is the dimensionless Morse potential. £j is a random number chosen from the range [0.5 — 1.5] 
and is modeling the random well depth of the interaction potential. The left lattice end is 
fixed and the right end has an arbitrary coordinate L > 0. The specific lattice elongation is 
e = L/(N + l). 

The true resonance condition i^ k m k u k = 0; {m} G Z Q), can be achieved at some elonga- 
tion values (apparently dependent on the set {&} of random values in (CQ)) as it is demonstrated 
in Fig.[TJ The most long- wave mode (k = 1) with the energy E$ = 10 -5 was initially excited and 
the energies of vibrational modes under adiabatically slow elongation are shown. Small number 
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FIG. 1: Evolution of the vibrational modes energies vs. the specific elongation e. The most long-wave 
mode (k = 1) is initially excited with the energy Eq = 10~ 5 . Digits numerate the modes numbers. 
N = 10. Modes not shown (k = 5, 6, 9, 10) are not excited. 

of particles N = 10 and low excitation energy allow to avoid the resonances overlap. One can 
see that the energy exchange is observed at some values of e. For instance, pair resonances 
u)3 : uji = 3 (e ~ 0.101), u; 2 ■ u)\ = 2 (e ~ 0.154) are obvious. Moreover, triple resonances are also 
presented (see e.g. u 3 + co> 4 = u}$ at e ~ 0.255). The energy exchange rates and degrees depend 
on E Q and the elongation rate. The energy decrease in Fig. [T] is explained by the "softening" 
of the potential U in under the elongation [^j], i.e. under adiabatic deformations the mode 
energy Ei decreases as its frequency: Ei ~ u>i, and w ~ ^/g, where g - lattice rigidity. 

The resonance u>2 :oo\ = 2 at e ~ 0.154 is chosen for detailed investigation. But firstly we 
consider the two-oscillatory model, where the resonance behavior can be analytically described. 



B. Two-oscillatory model. 

We consider two nonlinearly coupled harmonic oscillators with the following hamiltonian 
(here the main results are briefly demonstrated, whereas the full analysis will be presented 
elsewhere): 

H = ^ {pI + Pi) + 2 ( K i <Zi 2 + ^ <7 2 2 ) + Ag? <? 2 , (2) 

where Ki = 1, k-i = 4, A is the coupling parameter, and the frequencies ratio is 1:2. 

After the variables replacements qi i— > — /uJi cos (fu and (2/i-/ 2 )/5 i-> /, (2/ 1 -/ 2 )/5 i-> 
J, (20i — 02 ) i-> @, (0i + 20 2 ) i-> $ and the averaging 10, [lj]] over the period 2-n/u 2 Eq. (j2]) 
transforms to 

#(/, J) = 57-7(7 + 2 J) V27- J cosG, (3) 
where 7 = A. 7 is the integral of motion and the time evolution is determined by the 



second term. The energy exchange is governed by the following hamiltonian 

H-5I 



H = — = (l + 2p)^2-pcosO, (4) 

7i 

where p = J /I. The initial conditions po = 2 and po = —1/2 correspond to the excitation of 
the oscillators iO\ — 1 and cj 2 = 2, correspondingly. The Hamiltonian equations for TC(p, 0) 
have the form: 

p = (1 + 2p) V2 -p sinG , 

7-6p (5) 
6 = — . cos B . 
2 v / 2^ r p 

1. Excitation of the "stable" oscillator. The initial conditions p(0) = —1/2, O(0) = 
correspond to the excitation of the oscillator u 2 = 2. In this case the solution of Eq. (j^J) is 

p(t) = -1/2, 

. exp(2 v / T0t) - 1 (6) 

sin O(i) 



exp(2V10t) + 1 

and p(t) = const, what means that the oscillator u> 2 has the infinite life-time and the oscillator 
Ui is not excited at all. 

2. Excitation of the "unstable" oscillator. The situation drastically changes if the os- 
cillator uj\ = 1 is initially excited. Since (p = 2, O = 0) is a singular point of (jSJ), the initial 
condition p(£ = 0) = po — 2 — <5 (<5 <C 1) should be chosen with the final limiting transition 

5 — > +0. Taking into account that © has the integral of motion cosO = —, — ^ — ^ ° 
S ... (l + 2p)V2^ 

one can get the following solution in the vicinity of p(t) ~ po : 



and pit) varies through a finite range in a finite time ~ y4 Ap/a(po) and the phase trajectory 
goes away from the singular point po — 2 over a finite distance. Then, according to (JSJ), p(i) 
decreases approaching the value p = — 1/2. It means that the initially excited oscillator oj\ loses 
its energy in a finite time and completely gives the energy up to the higher frequency oscillator. 
This allows to specify the oscillator uj\ as an unstable, and the oscillator uj 2 as a stable one. 
3. "Mixed" initial conditions. Below we briefly consider the "mixed" initial condition 
(p = — 1/2 + v, 0) when the stable oscillator u 2 is initially excited with small [y <C 1) "weight" 
of the low-frequency oscillator. In this case the energy exchange is also observed and the 
life-times of both oscillators in excited states can be estimated as: 



t > J_ ln ( v^+l 

t 2 > —= In ■ 



(8) 



l + 2p 

and it follows that if the stable oscillator uo 2 is initially excited with a small addition v <C 1 of 
the oscillator u\, then the energy exchange is observed and the life-time of the unstable lower- 
frequency oscillator t\ is finite and does not depend on u, whereas the life-time t 2 of the stable 
higher-frequency oscillator increases ~ — ln(z/) as v — > +0. 



4. Resonance detuning. Obviously, the considered exact resonance ^2 = 2^ rarely occurs 
and the situation of resonance detuning is of particular interest. Let K\ — (1 — ft) and k 2 = 4 
in Eq. (J2J), where /i < 1 means a small resonance detuning. The case of the initial excitation 
of the stable oscillator is trivial: the equations of motion of §2§ are: 

q\ = -(1 - ft) qx - 2A qiq 2 ] q 2 = -4 q 2 - A q\ (9) 

and if qi(t = 0) = 0, q 2 (t = 0) ^ then the first oscillator is not excited at all, and the 
initially excited second oscillator lives infinitely long independently on the vicinity to the exact 
resonance. 

The case of the initial excitation of the unstable oscillator u)\ requires more accurate analysis. 
It is convenient to make the variables replacements: (q\ cos a — q 2 sin a) i— > y\ and (gi sin a + 
q^cosa) I— > y 2 . Then by neglecting the terms of the order of A/x 1 one can choose a such 
that the true resonance problem is fulfilled in new variables {yi,y 2 }. But this problem was 
considered above. And finally the life-times of both oscillators are estimated as: 

T X (fi) = — L= In [ ) ; 

vTO \V3-l) (10) 

Thus if the unstable oscillator u\ is initially excited, then its life-time is a finite value and 
doesn't depend on the closeness to the exact resonance point, while the life-time of the stable 
oscillator increases as ~ — ln(/x) at fi — > +0. 

Below we numerically compare the kinetics and the energy exchange degree for the Morse 
lattice and the two-oscillatory model. 



III. NUMERICAL SIMULATIONS OF RESONANCE INTERACTION IN THE 
MORSE LATTICE AND TWO-OSCILLATORY MODEL. 

A. Comparison of the Morse lattice and two-oscillators model. 

Our primary goal is to investigate the rate and the energy exchange degree near the resonance 
for two most long-wave modes (k = 1, 2) in the Morse lattice depending on the model parameters 
(resonance detuning, excitation energy). The comparison is made with the two-oscillatory 
model and one-to-one correspondence between vibrational modes k — 1, 2 in the Morse lattice 
and lower and higher frequency oscillators is assumed. 

Note, that because of the nonlinearity of both models, the mode/oscillator frequencies de- 
pend on the excitation energy Eq. And the true resonance conditions are fulfilled for some 
values Eo(Eq) and fio(E ) ^ 0. It is convenient to represent the actual values of e and fi as 
e = e (E ) +Si and ft = fto(E ) +ft\, where £\ and fti are the measures of resonances detuning. 

We compare the resonance 2 u\ = uj 2 at e ~ 0.154 (see Fig. [1]) in the Morse lattice and the 
resonance in the two-oscillatory model with K\ — 1 — ft, K 2 — 4; parameter A for definiteness 



is chosen A = 0.3. The exact resonance values /xq and Eq are determined as values, where the 
largest life-times of the stable mode/oscillator are achieved in numerical simulation. These 
values are e {E ) = 0.153 808 2813623 and fj. {E ) = 3.374 961 52 • 10" 8 , and E = 10^ 6 in both 
cases. All digits in Eq(E ) are accurate: variation of the last (13th) digit noticeably diminishes 
the vibrational life-time. 

Thus one can analyze and compare the interaction of vibrational modes and oscillators vs. 
detuning parameters e\ and fx\ in both systems. 

B. Numerical results for the Morse lattice and oscillatory model. 

The temporal behavior of modes/oscillators energies is shown in Fig. [2] at different values 
of detuning parameters E\ and fi\. Here the unstable mode/oscillator are initially excited with 
Eq = 10~ 6 . The general scenario of the kinetics and the energy exchange degree is very similar 
for both systems. 

There is no interaction between the vibrational modes/oscillators far away from the res- 
onance conditions (|ei| > 10~ 2 and \fii\ > 10~ 3 ). There exists partial and symmetrical 
energy exchange in the middle range of the detuning parameters 10~ 4 < < 10~ 2 and 
10 -6 l/^i I ~ 10~ 3 (Figs. Eh-b). As the resonance point is approached (by diminishing |ei| and 
the energy exchange degree and the modes asymmetry increases. And at a certain detun- 
ing parameters values the full energy exchange is achieved (Fig. [2b). Further approaching to the 
resonance point results in the increase of the stable mode/oscillator life-times, whereas the life- 
times of the unstable mode/oscillator do not change (Figs. [2H-e). And finally the life-times of 
the stable mode/oscillator achieve their maximal values at the resonance point |ei| = \fii\ = 0.0 
(Fig. ISO- 

The dependence of the stable mode/oscillator life-times vs. values £i//ii is also the same 
for both models. The dependence is logarithmic (see Fig. [3h), as predicted by ( |T0|) . We also 
found that the unstable oscillator life-time T\ with very good accuracy varies as T\ ~ 1/A in a 
wide range (10~ 2 < A < 10 2 ), where A is the coupling parameter in (j2J). The mode/oscillator 
life-times are defined as a time intervals for full energy exchange. For example, T2 is a time 
interval between points A and B in Fig. EH. Similarly T\ is a time interval between points B 
and C in Fig. EH. 

Fig. [3b demonstrates the energy exchange degree AE/E vs. £i//Ji for both models in the 
middle range of the resonance detuning parameters (10~ 4 < < 10~ 2 and 10 -6 < < 
10 -3 ). One can see the extreme sensitivity of life-times and energy exchange degree to the 
variations of detuning parameters. 

As pointed above, the exact values for the onset of resonances depend on the excitation 
energy Eq. Hence the energy exchange kinetics should be also sensitive to Eq. Fig. H] shows 
the dependencies of the life-times of both modes in the Morse lattice and both oscillators 
vs. the energy of initially excited unstable mode/oscillator E . The extreme sensitivity of 
stable mode/oscillatory life-times in the vicinity of exact resonance (eq = 0.153 808 2813623 
and fi Q = 3.374 961 52 ■ 10~ 8 at E = 10" 6 ) is observed. 
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FIG. 2: Temporal modes/oscillator energies evolution at different values of the resonance detunings 
for the Morse lattice (blue color) and the two oscillators model (red color). Stable mode/oscillator 
- solid lines, unstable mode/oscillator - dashed lines, a) e\ = —2.380 828136 23 • 10~ 2 , fi\ = 10 -3 ; 
b) e x = -8.082 813 623 • 10" 4 , m = 4 ■ 10~ 5 ; c) ei = -5.828136 23 • 10~ 5 , m = 1.9 • 10~ 6 ; d) 
E\ = —8.136 23 • 10~ 8 , fj,i = 10 -9 ; e) e\ = 0.0, fii = 0.0. In all cases the unstable mode/oscillator are 
initially excited with Eq = 10~ 6 . Note the enlarged time scale in panel e). 

Thus, in the case of the initial excitation of unstable mode/oscillator the resonance behavior 
is the same in both models and is perfectly described by the analytical results in Sec. Ill Bl But 
the models have different features in the case of stable mode/oscillator initial excitation. In 
the oscillatory model the initially excited stable oscillator lives infinitely long independently 
on the closeness to the resonance point. The stable mode in the Morse lattice has finite life- 
time and the energy exchange is observed, the resonance behavior being independent on the 
resonance detuning. This different features of the models is caused by the addition terms in 
the Morse potential expansion which influence the energy exchange. Actually, an addition of 




FIG. 3: (a) The dependence of the stable mode/oscillator life-times T2 on the detuning parameter 
values for the Morse lattice and the oscillatory model in the vicinity of the exact resonance point; 
(b) the energy exchange degree vs. the detuning parameter for the Morse lattice and the oscillatory 
model. The unstable mode/oscillator are initially excited with Eq = 10~ 6 . Blue color - Morse lattice, 
red color - oscillatory model in both panels. 
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FIG. 4: The dependence of modes/oscillators life-times Tj on the energy of the initially excited 
unstable mode/oscillator Eq. Values £0 = 0.153 808 2813623 and /j,q = 3.374 96152 • 10~ 8 correspond 
to the exact resonance point at Eq = 10 -6 (ei = fix = 0.0). Morse lattice - blue color, oscillatory 
model - red color; solid lines - stable mode/oscillator, dashed lines - unstable mode/oscillator. 

the term A qi q\ in (j2J) results in the mode energy exchange and the finite life-time of the stable 
vibrational mode. 

The case of "mixed" initial conditions when the stable mode/oscillator is excited with a small 
addition of the unstable mode/oscillator is also analyzed and results are in good agreement with 
analytical predictions (see ( jTUl) ). Fig. [5] shows the dependence of the stable mode/oscillator 
life-times vs. the contribution of the unstable mode/oscillator energy into the total initial 
energy Eq. Here the initial energies of the stable and unstable modes/oscillators are (1 — v)Eq 
and vEq, correspondingly. 

Note that the observed resonance interaction of the vibrational modes is the general phe- 
nomenon in random lattices under the deformations at low temperatures. Actually an expan- 
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FIG. 5: The dependence of the stable mode/oscillator life-times T2 on the contribution of the unstable 
mode/oscillator energy v into the total initial energy Eq = 10 -6 . The initial energies of the stable 
and unstable modes are (1 — u)Eq and vEq, correspondingly; e\ = fi± = 0.0. Blue color - T2 for the 
Morse lattice, red color - T2 for the oscillatory model. 

sions of any interaction potentials have the form U ~ x 2 ± 71 3 . The minus sign means that 
the potential softens under elongation (as in the Morse lattice). On contrary the plus sign (as 
in the a-FPUT lattice) means that the potential softens under contraction. This is the reason 
which allowed to observe the similar resonance behavior in the random FPUT lattice under 
contraction. 



In conclusion, we briefly summarize the main results. 

1. The vibrational modes resonance was found in the random Morse lattice. Two features allow 
to observe the resonances: the lattice randomness (random values of the interaction potential; 
the random masses also can be used), and the lattice deformation. At some values of the specific 
elongation the exact resonance condition m>i Ui = are be fulfilled. The resonance 2u\ = 002 
of two most long-wave vibrational modes was chosen for the detailed analysis. Small number 
of particles (N = 10) and low excitation energy allowed to avoid the resonances overlap. 

2. Simple two-oscillators model with inharmonic coupling was also considered. The analytical 
description of the kinetics and modes energy exchange vs. the model parameters (resonance 
detuning, specific energy) was performed. 

3. Surprisingly good agreement was observed for the kinetics and energy exchange rate for the 
Morse lattice and simple model of two nonlinear coupled harmonic oscillators. 

4. Modes/oscillatory life-times and the energy exchange degree are extremely sensitive to the 
resonance detuning parameters. 

5. The resonance behavior is a general feature of random nonlinear lattices at deformations 



IV. CONCLUSIONS. 



and it was also observed in the random a-FPUT lattice. 



The observed resonance phenomenon can not throw some light on the problem of energy 
equipartition and onset of the chaos because of small values of N and low specific energy, but it 
can serve as a starting point for the analysis of initial stages of vibrational modes interaction. 
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